3.4.1 Ts sources (MGRsD, MGRsD_TAIR, From paper, Rs_Ts_relationship)

First, we tested the effect of different soil temperature sources on the Rs_annual_bahn vs Rs_annual relationship:

• From MGRsD means mean annual soil temperature (amst) are from a global monthly soil respiration database, each site has more than 12 months measured soil temperature read from original papers.
• From paper means amst were reported from the original paper (table, figures, or description).
• Partly from TAIR means: some studies did not measure soil temperature all year, for those months, we predict soil temperature based on monthly air temperature (Tsoil = 2.918 + 0.829*Tair). This model was developed based on the sites which have >= 12 months soil temperature measurements.

There are 67 records for which I cannot get the soil temperature information through the above three methods. In these cases, based on the Rs_Ts_relationship and reported Rs_annual, I calculated the amst.

• The calculated amst was then compared with the annual Tair, if they are well matched (error < 5%), calculated mast were used.
• When calculated amst and annual Tair do not match, this usually indicates a potential problem, and then I went back to the manuscript and checked.
• Whenever a paper reported annual mean Ts, I compared the reported mast and estimated mast based on the Rs_Ts_relationship, and found they are well matched.

Generally, Ts sources do not have clear effects on the Rs_annual_bahn and Rs_annual relationship.

3.4.2 Annual Rs or Ts coverage effect

Second, we tested whether Ts and Rs coverage (e.g., 0-0.5 means Rs or Ts only measured less than 6 months, versus the entire year), but found that Ts and Rs coverage do not have significant effects on the Rs_annual_bahn vs Rs_annual relationship.

3.4.3 Effect of maximum allowed divergence between global climate data set and site-specific air temperature

Third, we tested whether maximum allowed divergence between global climate data set and site-specific air temperature affect the Rs_annual_bahn vs Rs_annual relationship.

As we throw out data points with high divergence, R2 and RSE go up and down inconsistently and by small amounts, suggesting that the Tair divergence do not have a consistent effect.

3.4.4 Effect of maximum allowed divergence between annual precipitation from paper and Del

Fourth, we also tested whether the maximum allowed divergence between the global climate data set and site-specific precipitation affects the Rs_annual_bahn vs Rs_annual relationship. In other words, does a bias in the global data affect things?

As we throw out data points with high divergence, R2 and RSE showed no large or consistent changes, suggesting that the precipitation divergence does not have a large effect.

We also compared: * MAT from U. Delaware (MAT_Del) and MAT reported from the papers (a) * TAnnual from U. Delaware (TAnnual_Del) and annual temperature from papers (study_temp) (b) * MAP from U. Delaware (MAP_Del) and MAP reported from the papers (c) * PAnnual from U. Delaware (PAnnual_Del) and annual precipitation from papers (study_precip) (d)

In general, the temperature and precipitation from the University of Delaware climate data matched the data reported from publications well. This supports the idea that the divergence between global climate data set and site-specific precipitation/temperature has little to no effect.

We tested the effect of precipitation and temperature variability (quantified by standard deviation from 1961 to 2014), using multiple linear regression, with divergence as catergorical indicator. [[JJTODO: this needs to be clearer. sd of what exactly? Annual means and sums?]]

We found that they have no significant effect on the Rs_annual_bahn vs Rs_annual relationship.

3.4.5 Test ecosystem type

The Rs_annual_ban vs Rs_annual relationship varies among different ecosystems.

• For example, agriculture has lower slope but wetland has higher slope.
• However, it is unlikely that the data from agriculture and wetland shift the overall regression between Rs_annual_bahn vs Rs_annual away from 1:1 line.
• When we remove the Ag data, the Rs_annual_bahn vs Rs_annual regression still differs from 1:1 line.

• When we plot Ag alone, we see the Rs_annual_bahn vs Rs_annual regression in Ag does not greatly differ from the rest.
• We also examined the influential outlier points (cooks.distance > 0.5). When the outliers were removed, the regression showed no difference, indicating that outliers do not have large effects.
• Similar conclusion for wetland data.

3.4.6 Test Rs measurement method

Different measurement methods do not affect the Rs_annual_bahn vs Rs_annual relationship.

3.4.9 Drought effect?

In their paper, Bahn et al. (2010) reported that drought stress significantly affected the relationship. Using our new datasets we found that:

• MAP effects on the Rs_annual_bahn vs Rs_annual relationship is very limited (p ~~ 0.05).
• The slope and intercept changes do not followed a clear pattern as MAP changes.
• MAP is not a good drought index.

We then tested standardized drought index (SPI).

• SPI significantly affects Rs_annual_bahn vs Rs_annual relationship; as SPI decreases, the slope decreases, means the Bahn approach tends to overestimate Rs_annual under drought condition.

Since SPI compares annual precipitation at a site with average precipitation over a period (we used 1964-2014), it can not describe spatial drought in drought.

We thus also used another drought index, the Palmer Drought Severity Index (PDSI), to characterize spatial effects.

• PDSI significantly affects the Rs_annual_bahn vs Rs_annual relationship: as it becomes drier, the slope tend to decrease, meaning that the Bahn method overpredicts the annual flux. This is consistent with the SPI finding above.

4.2 Predicting Rs_annual from air temperatures

Since high resolution soil temperature is still lacking, and/or has lower accuracy than air temperature data, we want to test whether we can use Rs at annual mean air temperature (amat) or mean annual temperature (mat) to predict Rs_annual. [JJTODO: clarify the difference between these]

• The regression between Rs_annual_bahn_amat and Rs_annual falls away from the 1:1 line, with an intercept significantly different from 0 (p<0.001) and slope significantly different from 1 (p<0.001).
• Note that we tested both new1 and new2 models.
• Using amat or mat show no big difference, but it is surprising that using mat is slightly better than using T_Annual. [JJTODO: be consistent in terminology]

We detected 2 outliers in the Rs_annual_bahn vs Rs_annual regression. * Remove these two outliers significantly improved the model (slope changed from 0.78 to 0.87, intercept decreased from 222 to 157, however, p values for slope and intercept are still < 0.001).

## Warning: Removed 2 rows containing non-finite values (stat_smooth).

## Warning: Removed 2 rows containing missing values (geom_point).

## Warning: Removed 1 rows containing non-finite values (stat_smooth).

## Warning: Removed 1 rows containing missing values (geom_point).

[JJTODO: this is very unclear] * We adjusted the air temperature by TAnnual_adj = 2.918+0.829TAnnual MAT_adj = 2.918+0.829MAT The results are better, but still not resolve the problem (p<0.05)

4.3 Re-simulate a model (new3)

We re-calculated soil respiration at annual mean air temperature (Rs_amat, i.e., using air temperature rather than soil temperature to calculate soil respiration).

• Then we re-simulated the relationship between Rs_annual and Rs_amat.
• Rs_annual = 729.09225 * (Rs_amat ^ 0.46535) + 89.7789 * spi – Medeterrean
• Rs_annual = 588.618 * ( Rs_amat ^ 0.65022 ) + 22.59026 * spi + 11.29775*pdsi – exclude Mediterranean sites
• If we update the model, Rs_annual_bahn_amat can represent Rs_annual.